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Abstract  Nitric  oxide  planar  laser-induced  fluorescence 
(NO  PLIF)  imaging  is  demonstrated  at  a  10- kHz  repetition 
rate  in  the  Calspan  University  at  Buffalo  Research  Center’s 
(CUBRC)  48-inch  Mach  9  hypervelocity  shock  tunnel  using 
a  pulse  burst  laser-based  high  frame  rate  imaging  system. 
Sequences  of  up  to  ten  images  are  obtained  internal  to  a 
supersonic  combustor  model,  located  within  the  shock  tun¬ 
nel,  during  a  single  ~  10-millisecond  duration  run  of  the 
ground  test  facility.  Comparison  with  a  CFD  simulation 
shows  good  overall  qualitative  agreement  in  the  jet  pene¬ 
tration  and  spreading  observed  with  an  average  of  forty 
individual  PLIF  images  obtained  during  several  facility  runs. 

1  Introduction 

The  need  for  optical  diagnostic  methods  suitable  for 
application  in  high-enthalpy  hypervelocity  short-duration 
facilities  is  especially  acute  due  to  the  unique  capabilities  of 
these  facilities  for  large-scale  testing  of  hypersonic  pro¬ 
pulsion  technologies  at  duplicated  flight  conditions  as  well 
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as  the  fundamental  physics  of  hypervelocity  flows.  How¬ 
ever,  while  the  development  of  non-intrusive  diagnostic 
techniques  has  experienced  much  recent  progress,  several 
challenges  remain  due  to  the  inherent  constraint  of  limited 
optical  access  and  typically  short  (order  100’s  of  micro¬ 
seconds  to  tens  of  milliseconds)  steady  state  run  times. 
Limitations  on  the  number  of  facility  “shots”  per  day 
(typically  two  to  three  in  the  current  facility)  are  another 
consideration.  Despite  these  challenges,  several  studies 
incorporating  advanced  optical  diagnostic  methods  in  such 
facilities  have  been  recently  reported.  For  example,  nitric 
oxide  molecular  tagging  velocimetry  (NO  MTV)  (Danehy 
et  al.  2009),  Doppler-shift-based  velocimetry  (Danehy  et  al. 
2001;  Danehy  and  O’Byme  1999),  and  NO  planar  laser- 
induced  fluorescence  (PLIF)  measurements  of  number 
density  and  thermometry  (O’Byme  et  al.  2002)  have  been 
performed  in  the  Australia  National  University’s  T2  free- 
piston  shock  tunnel.  More  recendy,  two-component  Doppler- 
shift-based  velocimetry  has  been  demonstrated  in  the 
University  of  New  South  Wales’  T-ADFA  free-piston 
shock  tunnel  (Hruschka  et  al.  1996).  As  another  example,  a 
novel  fuel  plume  imaging  diagnostic,  utilizing  Mie  scat¬ 
tering  from  in  situ  generated  silica  particles  in  combination 
with  a  long  (100  microsecond)-pulsed  Alexandrite  laser, 
has  been  developed  and  applied  to  injector  studies  in  the 
NASA-HYPULSE  facility  at  General  Applied  Sciences 
Laboratory  (GASL)  (Tsai  et  al.  1996).  NO  number  density 
measurements  by  UV  absorption  and  NO  PLIF  have  also 
been  reported  in  the  HYPULSE  facility  (Roberts  et  al. 
1994).  Several  diagnostic  advances  have  also  been  reported 
in  studies  at  the  Calspan  University  at  Buffalo  Research 
Center  (CUBRC)  LENS  I  and  LENS  II  Hypersonic  Shock 
Tunnels,  including  measurement  of  temperature  and 
velocity  by  water  vapor-based  diode  laser  absorption  (Wehe 
et  al.  1988);  NO  number  density,  temperature  and  velocity 
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by  NO-based  diode  laser  absorption  (Parker  et  al.  2006, 
2007);  and  aero-optic  distortion  by  Phase  Shift  Interfer¬ 
ometry  utilizing  infrared  scene  generation  and  recording 
techniques  (Holden  et  al.  1993;  Holden  and  Parker  2002). 
Finally,  Barker  et  al.  (1998)  have  reported  a  novel  flow 
tagging  velocimetry  diagnostic  for  use  in  high-enthalpy 
short-duration  facilities  based  on  laser-enhanced  ionization 
of  atomic  sodium;  Allen  et  al.  (1993)  have  reported  single- 
laser  pulse  NO  PLIF-based  thermometry  in  a  high-enthalpy 
shock  tunnel  facility  located  at  PSI  corporation;  and 
McMillan  et  al.  (1993)  have  reported  two-line  NO  PLIF 
temperature  imaging  of  a  transverse  jet  in  a  supersonic 
cross-flow.  All  of  these  imaging  diagnostics  have  been 
performed  with  conventional  ~  10-Hz  repetition  rate  pulsed 
laser  sources,  which  limits  data  accumulation  to  one  image 
per  facility  run. 

More  generally,  since  its  inception  more  than  twenty-five 
years  ago  (Dyer  and  Crosley  1982;  Kychakoff  et  al.  1984), 
PLIF  imaging  has  been  developed  into  a  powerful  technique 
for  qualitative  and  quantitative  visualization  in  a  wide 
variety  of  flow  and  combustion  environments.  However,  the 
ability  to  capture  unsteady  flow  dynamics  at  high  speed  is 
severely  constrained  by  the  relatively  low  repetition  rates 
(order  10  Hz)  of  typical  PLIF  systems.  Such  instrumenta¬ 
tion  could  only  acquire  a  single  image  in  one  shock  tunnel 
run.  We  have  been  developing,  over  the  last  several  years, 
an  ultra-high  frame  rate  PLIF  imaging  system  (up  to  1  MHz 
frame  rate)  which  combines  a  “burst  mode”  Nd:YAG- 
based  tunable  optical  parametric  oscillator  (OPO)  system 
(Jiang  et  al.  2008;  Jiang  and  Lempert  2008)  with  a  variety 
of  high  frame  rate  cameras,  including  CCD,  ICCD  and 


intensified  CMOS  architectures.  In  addition  to  enabling 
capture  of  the  temporal  dynamics  of  hypervelocity  flows  at 
frame  rates  up  to  1  MHz,  by  operation  at  lower  frame  rates 
of  the  order  of  10-20  kHz,  the  burst  mode  imaging  system 
provides  the  unique  ability  to  obtain  sets  of  8-10  full  planar 
images  in  a  single  short-duration  facility  test.  The  work 
reported  here  was  motivated  by  the  desire  to  demonstrate 
such  potential.  Specifically,  we  demonstrate  the  ability  to 
obtain  ten  high  resolution  NO  PLIF  images  during  a  single 
run  of  the  CUBRC  48-inch  shock  tunnel,  operating  at 
nominal  Mach  9  free  stream  conditions.  Images  of  injector 
plume  mixing  are  obtained  in  a  large-scale  Multidisciplin¬ 
ary  University  Research  Initiative  (MURI)  high-speed 
combustion  duct  in  which  NO-seeded  helium  is  injected 
into  a  pure  N2-free  stream.  Comparison  with  a  CFD  simu¬ 
lation  shows  good  overall  agreement  in  fuel  penetration  and 
downstream  spreading  of  the  plume. 

2  Experimental 

2.1  Laser  system  description 

The  burst  mode  pump  laser,  illustrated  in  Fig.  1,  has  been 
described  in  detail  previously  (Jiang  et  al.  2009;  Thurow 
et  al.  2004)  and  will,  therefore,  be  only  summarized  here. 
Briefly,  a  continuous  wave  (cw)  diode-pumped  Nd:YAG 
ring  laser  serves  as  the  primary  oscillator,  the  output  of 
which  is  preamplified  in  a  double-pass  variable  pulse  width 
(from  0.3  to  2.0  ms)  flashlamp-pumped  pulsed  amplifier. 
The  amplifier  rod  is  10  cm  in  length  by  6.4  mm  in 


Fig.  1  Schematic  diagram  of 
five-amplifier  configuration 
pulse  burst  Nd:YAG  laser 
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diameter,  and  is  wedged  at  a  few  degrees  to  mitigate  self¬ 
lasing.  The  resulting  preamplified  output  is  formed  into  a 
“burst”  train  of  short  pulses  using  a  custom,  dual  Pockels 
cell  “sheer.”  More  specifically,  the  slicer  produces  burst 
sequences  of  8-10  pulses  at  1,064  nm,  each  with  duration 
of  —  10  ns,  and  pulse  energy  of  ~  10-9  J  (1  nJ).  The  series 
of  pulses  is  further  amplified  in  a  series  of  four  additional 
flashlamp-pumped  amplifiers,  with  6.4,  6.4,  9,  and  12.7 
mm  diameter  x  10-cm  long  rods,  respectively.  This  con¬ 
figuration  is  somewhat  different  from  the  seven  amplifier 
configuration  which  has  been  used  previously  for  operation 
at  burst  frequencies  of  up  to  1  MHz  (Jiang  et  al.  2011).  As 
will  be  discussed  below,  the  desire  to  obtain  the  highest 
possible  spatial  resolution  dictated  the  use  of  a  10-kHz 
frame  rate-intensified  CMOS  camera,  rather  than  the 
1  MHz  framing  CCD  camera  used  in  most  previous  high¬ 
speed  studies.  The  lower  frame  rate  in  turn  dictated  a 
correspondingly  low  (10  kHz)  burst  rate  and  corresponding 
burst  envelope  of  1  ms  (10  pulses  at  100  microsecond 
spacing).  Due  to  very  significant  gain  refilling  between  the 
individual  pulses  in  the  burst,  only  five  amplifiers  are 
required  to  reach  the  high  energies  required  for  these 
measurements.  The  remaining  amplifier  stages  increase  the 
individual  average  pulse  energy  from  1  nanoJoule  to 
—  400  mJ  per  pulse  (4  x  108  gain)  at  1,064  nm. 

Note  that  amplifiers  cannot  simply  be  added  arbitrarily, 
as  system  performance  is  ultimately  limited  by  the  growth 
of  amplified  spontaneous  emission  (ASE)  in  the  forward 
direction.  Inclusion  of  a  stimulated  Brillouin  scattering 
(SBS)  phase  conjugate  mirror  (PCM)  is  a  well-known 
approach  for  improving  the  spectral  characteristics  of  a 
pulsed  laser  beam  (Ni  and  Kung  1996).  For  this  work,  the 
purpose  of  the  PCM  is  to  both  improve  the  on/off  contrast 
of  the  pulse  slicer  and  to  limit  ASE  in  the  forward  direc¬ 
tion.  Fluorinert  FC-75,  which  has  been  used  extensively  as 
an  SBS  fluid  (Jiang  et  al.  2008)  and  has  previously  per¬ 
formed  well  when  operating  the  burst  mode  laser  at  less 
than  ~  250  kHz  (Jiang  et  al.  2008),  was  used  for  this  work. 

The  fundamental  output  at  1,064  nm  is  converted  to 
second-  (532  nm)  and  third  (355  nm)-harmonic  wave¬ 
lengths  using  a  pair  of  KD*P  crystals.  The  third  harmonic 
is  then  used  as  the  pump  for  a  home-built  injection-seeded 
OPO  system.  The  OPO  output  is  sum-frequency  mixed 
with  residual  355  nm  from  the  pump  laser,  creating  high- 
frequency  bursts  of  tunable  output  in  the  vicinity  of 
226  nm,  which  are  used  for  NO  PLIF  imaging.  The  repe¬ 
tition  rate  of  the  burst  sequence  can  be  as  high  as  10  Hz, 
but  is  typically  limited  to  1—4  Hz  when  operating  with 
burst  envelopes  exceeding  —  0.5  ms,  due  to  thermal  load¬ 
ing  of  the  amplifier  rods.  Figure  2  shows  a  typical  10-pulse 
oscilloscope  trace  obtained  at  10-kHz  burst  rate  at  532  nm. 
The  average  individual  pulse  energy  is  —  150  mJ.  Note  that 
the  trace  was  obtained  with  a  relatively  low  sampling  rate 


Time  (|xs) 

Fig.  2  Typical  ten-pulse  laser  burst  trace  at  532  nm  with  100-micro- 
second  inter-pulse  spacing  (10-kHz  repetition  rate).  Average  pulse 
energy  is  ~  150  mJ 

oscilloscope  that  required  the  detector  output  to  be  delib¬ 
erately  broadened  in  time  (using  a  load  resistor  and  the 
intrinsic  cable  capacitance). 

For  spectrally  narrowed  imaging,  the  OPO  was  operated 
with  injection  seeding  (Jiang  et  al.  2008)  using  an  external 
cavity  diode  laser  (ECDL)  manufactured  by  Sacher,  Inc. 
The  OPO  “idler”  beam  is  seeded,  which  results  in  line- 
narrowed  output  for  both  the  idler  and  signal  waves,  with 
time-averaged  linewidth  of  —  300  MHz,  based  on  mea¬ 
surements  using  a  commercial  Fizeau  wavemeter.  For  the 
seeded  measurements  performed  in  this  work,  the  seed  laser 
was  tuned  to  822.044  nm  (vacuum  wavelength),  measured 
with  a  HighFinesse  WS-6  Fizeau  interferometer  wavemeter. 
Mixing  the  output  of  the  OPO  with  the  third  harmonic  of  the 
Nd:  YAG  laser  produced  UV  output  at  226.244  nm,  exciting 
the  Qi  |(5.5)  +  QP21  (5.5)  transition  of  the  NO  A-X  (0,0) 
band  (using  the  notation  of  Reisel  et  al.  (1992)).  The 
spectral  linewidth  at  226  nm,  estimated  from  measurements 
of  the  linewidth  of  the  822  nm  and  355  nm  output, 
is  —  0.5  GHz,  which  is  approximately  equal  to  the  NO 
pressure-broadened  linewidth  for  static  pressures  of 
—  0.5  bar  and  static  temperatures  of  —400  K,  which  are 
typical  values  for  the  central  portion  of  the  imaged  region  of 
the  flow  (See  Figs.  8,  9).  This  transition  was  selected  based 
on  the  fact  that  the  nominal  static  temperature  of  the  fuel 
plume  was  predicted  to  vary  between  —  1 00  K  at  the  nozzle 
exit  (i.e.,  the  near  wall  region  of  the  jet)  and  —400  K  after 
mixing  with  the  free  stream  flow.  Figure  3  shows  a  set  of 
226-nm,  10  KHz  burst  traces.  The  average  individual  pulse 
energy  is  —0.2  mJ/pulse. 

The  UV  beam  had  a  circular  cross-section  with  a  5  mm 
diameter,  and  it  was  directed  upwards,  using  a  series  of 
mirrors,  to  the  top  of  the  facility.  After  propagation 
(  — 5  m)  to  the  shock  tunnel,  however,  the  beam  had 
diverged  to  —  10  mm  in  diameter.  The  laser  beam  was  then 
formed  into  a  thin  sheet  using  the  combination  of  a 
cylindrical  lens  and  a  spherical  lens.  A  positive  300-mm 
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Fig.  3  Typical  ten-pulse  laser  burst  trace  at  226  nm  with  100-micro- 
second  inter-pulse  spacing.  Average  pulse  energy  is  ~  0.2  mJ 

focal  length  cylindrical  lens  focused  one  axis  of  the  beam 
which  then  propagated  600  mm  and  passed  through  a 
1 -meter  focal  length,  50-mm  diameter  spherical  lens.  This 
pair  of  lenses  increased  the  total  width  of  the  beam 
to  ~  100  mm.  Of  this  100  mm,  the  central  50  mm  was 
used  for  the  actual  imaging.  Using  only  the  central  portion 
resulted  in  improved  sheet  uniformity  and  also  provided 
some  accommodation  to  tunnel  motion,  which  was  known 
to  translate  parallel  to  the  primary  flow  axis  during  facility 
operation.  Since  the  beam  delivery  optics  was  not  fixed 
with  respect  to  the  tunnel,  this  motion  could,  in  principle, 
result  in  clipping  of  the  sheet  at  the  facility  entrance  win¬ 
dow.  In  practice,  this  was  not  found  to  be  a  significant 
issue.  The  final  UV  laser  sheet  was  ~0.6  mm  thick. 

2.2  48-inch  shock  tunnel  facility  and  pulse  burst  laser 
facility  integration 

The  48-inch  shock  tunnel  facility  is  the  smallest  of  three 
such  hypersonic  facilities  located  at  CUBRC  in  Buffalo, 
N.Y.,  which  can  provide  duplicated  velocity/altitude  con¬ 
ditions  for  combustion  studies.  Fundamental  hypervelocity 
studies  from  high-enthalpy  reentry  physics  to  interceptor 
flow  dynamics  and  turbulent  combustion  chemistry  in  full- 
scale  combustion  ducts  have  been  studied  with  these 
facilities.  The  48-inch  shock  tunnel  can  generate  total 
enthalpy  conditions  in  excess  of  4  MJ/kg  required  for 
realistic  temperature  combustion  studies  at  Mach  8.  The 
steady  state  test  times  of  5-10  ms  at  this  high-enthalpy 
condition  provide  sufficient  flow  lengths  to  stabilize  com¬ 
plex  interaction  regions  and  the  combusting  flows  behind 
angled  wall  jets  in  a  flight  scale  (HIFiRE  II)  combustor. 
With  different  contoured  nozzles  and  throats,  the  48-inch 
shock  tunnel  can  support  test  flows  in  the  Mach  number 
range  from  6  to  20.  For  the  proof  of  concept  imaging 
experiments  described  in  this  work,  a  nominal  Mach  9 
nozzle/throat  configuration  was  used  to  produce  the  pure 
nitrogen  inlet  test  gas. 

The  MURI  combustion  duct  is  shown  supported  in  the 
test  section  of  the  48-inch  shock  tunnel  in  Fig.  5.  This 


combustion  duct  has  dimensions  of  40  cm  (height)  x 
40  cm  (span)  x  200  cm  (length).  The  rectangular  com¬ 
bustion  section  within  this  model  is  approximately  10  cm 
(span)  x  3  cm  (height)  x  30  cm  in  length  where  the  Mach 
number  is  approximately  3  based  on  CFD  computations. 
Five  0.254-cm  diameter  sonic  fuel  injectors  are  located  on 
the  bottom  wall,  spaced  6  mm  apart  and  oriented  at  30° 
with  respect  to  the  model  floor.  All  five  injectors  are  fueled. 
These  injectors  operate  at  an  injection  pressure  of  nomi¬ 
nally  500  psi.  For  all  images  to  be  shown  below,  the 
injected  gas  is  5  %  nitric  oxide  (NO)  in  helium,  which  is 
selected  as  a  non-combustible  surrogate  for  hydrogen,  in 
order  to  study  mixing. 

The  laser  optical  system  is  mounted  on  a  portable  optical 
table  so  that  it  could  be  employed  in  a  number  of  different 
wind  tunnel  facilities.  A  separate  portable  cart  was  assem¬ 
bled  for  the  flashlamp  power  supplies  and  other  electrical 
driver  components.  The  power  supplies  were  individually 
fused  to  two  standard  three-phase  220- V,  30-amp  breaker 
boxes.  As  illustrated  schematically  in  Fig.  4,  the  UV  output 
was  directed  through  two  fused  silica  windows,  located  at 
the  top  of  the  facility  and  test  section,  respectively.  The 
downward  propagating  laser  sheet  was  focused  directly  on 
the  bottom  surface  of  the  combustor  (there  was  no  exit 
window).  A  10  kHz  intensifier  (LaVision  HS-IRO)  and 
CMOS  camera  (Vision  Research,  Phantom  710)  (as  shown 
in  Fig.  5)  were  placed  close  ( —  150  mm)  to  the  test  section 
imaging  window,  resulting  in  the  object  plane  being 
approximately  0.75  m  from  the  camera  lens.  The  total 
camera  pixel  format  was  800  x  1,200.  As  alluded  to  earlier, 
this  camera  was  utilized  in  order  to  maximize  obtainable 
spatial  resolution  while  simultaneously  maintaining  an 
acceptable  field  of  view.  This  camera  was  used  in  combi¬ 
nation  with  a  singlet  plano-convex  UV  lens  (200-mm 
nominal  focal  length,  50  mm  diameter)  for  imaging.  The 
resulting  nominal  object  plane  spatial  scale  was  ~75 
microns/pixel,  with  actual  spatial  resolution  limited  in 
practice  by  the  LaVision  intensifier  and  the  performance  of 
the  plano-convex  lens.  This  imaging  was  accomplished 
using  a  nominal  1/3.4  magnification,  resulting  in  a  total 
object  plane  field  of  view  of  ~~  60  mm  (vertical)  x  85  mm 
(horizontal)  for  the  full  800  x  1,200  pixel  format  sensor. 
Note  that  the  intensifier  diameter  of  25  mm  and  the  solid 
collection  angle  of  —  f/15  were  sufficient  for  high  signal/ 
noise  instantaneous  NO  PLIF  imaging. 

Use  of  the  available  higher  frame  rate  ICCD  camera 
(Princeton  Scientific  Instruments,  PSI-IV),  which  has  a  pixel 
format  of  160  x  160,  would  have  resulted  in  a  spatial  res¬ 
olution  degraded  by  a  factor  of  approximately  ten.  Since  it 
was  initially  estimated  that  the  penetration  of  the  injector 
plume  into  the  free  stream  flow  would  in  some  cases  be  only 
on  the  order  of  1  -3  mm,  it  was  decided  to  optimize  maximum 
spatial  resolution  while  limiting  the  frame  rate  to  ~  10  kHz. 
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Fig.  4  Simplified  schematic 
diagram  illustrating  combustion 
duct  and  laser  propagation 
geometry.  Inlet  flow  is 
nominally  Mach  9  and  injection 
angle  is  approximately  30° 


Laser  Sheet 


UV  Window 


Fig.  5  Photograph  of  generic  scramjet  engine  model  and  intensified 
CMOS  camera  system 

During  the  shock  tunnel  experiments,  the  intensifier  gain 
was  fixed  at  ~  50  %  of  its  maximum,  and  the  exposure  gate 
set  to  ~500  ns.  A  custom-made  2-inch  diameter  trans¬ 
mission  filter  from  Layertec,  GmbH  was  placed  in  front  of 
the  plano-convex  lens  to  attenuate  the  elastic  laser  scatter¬ 
ing  while  transmitting  most  of  the  red-shifted  fluorescence 
(transmits  1  %  @  226  nm  and  >80  %  @235-280  nm). 
This,  in  combination  with  applying  black  paint  to  the  bot¬ 
tom  surface  of  the  combustor,  was  sufficient  to  eliminate 
stray  elastic  scattering  from  the  laser  sheet. 

Due  to  the  steady  state  test  time  of  —  5-10  ms,  a 
somewhat  complicated  scheme  was  required  in  order  to 
synchronize  the  tunnel  run  with  the  laser  and  camera  sys¬ 
tems.  In  particular,  the  CMOS  camera  system  required 
three  distinct  synchronization  signals,  and  the  pulse  burst 
laser  needed  to  be  continuously  triggered  at  1  Hz  during 
warm-up  prior  to  the  tunnel  run,  mostly  due  to  the  OPO  and 
sum-frequency  mixing  crystals  that  were  not  temperature 
stabilized.  The  timing  box  labeled  “SRS  Delay  Generator 
1”  (near  the  center  of  Fig.  6)  was  used  to  generate  this 
1  Hz  system  trigger.  It  triggered  a  second  SRS  delay 


Fig.  6  Block  diagram  illustrating  triggering  employed  for  shock 
tunnel  measurements 

generator,  which  in  combination  with  two  other  delay 
generators  provided  triggers  for  the  laser  flashlamp  power 
supplies  and  the  pulse  forming  pulse  slicer.  This  aspect  of 
the  synchronization/triggering  was  unchanged  from  our 
previous  work.  The  Phantom  CMOS  camera  required  a 
continuous  10  kHz  TTL  “clock”  pulse,  as  well  as  a  sepa¬ 
rate  10  kHz  burst  for  triggering  the  actual  image  acquisi¬ 
tion  series.  The  LaVision  intensifier  required  a  second 
10  kHz  burst  for  synchronization  to  the  CMOS.  The 
10  kHz  burst  output  of  the  primary  SRS  delay  generator 
was  inhibited  with  a  99.9  percent  duty  cycle  (1  ms  out  of 
each  1  s  cycle)  during  system  warm-up  (labeled  “Warm 
Up  Inhibit”).  Just  prior  to  a  tunnel  run,  this  inhibit  timing 
was  switched  to  a  second  timing  box  (labeled  “Shock 
Tunnel”)  which  was  triggered  from  a  pressure  transducer 
signal  indicating  that  the  primary  flow  had  been  estab¬ 
lished.  In  this  manner,  the  1-Hz  warm-up  was  maintained 
for  all  but  one  pulse  prior  to  the  facility  firing,  and  even  this 
was  manually  manipulated  to  give  a  near  1-Hz  interval 
between  the  warm-up  trigger  signals  and  that  of  the  facil¬ 
ity.  Note  that  this  procedure  resulted  in  a  maximum  of  ±50 
microsecond  jitter  between  the  triggering  of  the  camera 
system  and  the  start  of  acquisition,  which  is  negligibly 
small  compared  to  the  facility  run  time  of  5-10  ms. 
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3  Illustrative  imaging  results 

3.1  Prefacility  run  underexpanded  jet  images 

In  this  section,  we  present  a  sampling  of  illustrative  high- 
frequency  frame  rate  imaging  results  at  10  kHz.  Prior  to 
obtaining  images  from  shock  tunnel  runs,  sets  of  images 
were  obtained  from  the  model  in  the  identical  geometry  to 
that  which  would  be  used  for  tunnel  runs  but  with  the 
injector  operating  as  a  highly  underexpanded  jet  (i.e., 
emptying  into  a  nominal  vacuum  without  a  main  flow). 
This  provided  a  convenient  means  to  align  the  optical 
system,  determine  adequate  NO  seed  and  intensifier  gain 
levels,  test  the  triggering  and  warm-up  procedures,  and  also 
to  provide  an  approximate  assessment  of  the  image  quality. 
Figure  7  shows  a  typical  six  frame  sequence  of  highly 
underexpanded  jet  images  obtained  from  a  single  10  kHz 
laser  burst.  The  flow  is  NO/helium  supplied  through  the 
five  sonic  orifices  with  a  500  psi  plenum  pressure.  The  laser 
sheet  interrogates  the  flow  out  of  a  single  orifice.  The  field 
of  view  is  3.0  cm  x  6.5  cm,  with  corresponding  spatial 
scale  of  ~75  microns/pixel.  The  flow  is  propagating  from 
bottom  left  to  upper  right,  at  an  angle  of  ~  45°,  with  nozzle 
exit  plane  centered  at  the  (x,  y)  location  (1  cm,  —0.5  cm). 
While  admittedly  qualitative,  the  image  quality  is  quite 
reasonable,  illustrating  the  well-known,  and  rather  com¬ 
plex,  structure  of  the  underexpanded  jet.  In  particular,  the 
characteristic  “barrel”  shock,  with  corresponding  discon¬ 
tinuous  increase  in  density,  is  clearly  evident,  as  is  some 


irregular  (possibly  turbulent)  structure,  particularly  down¬ 
stream  of  the  Mach  disk  (normal  shock). 

3.2  Facility  run  images 

Figure  8  shows  an  example  eight-frame  image  sequence 
(top  left  to  bottom  right)  obtained  at  10-kHz  frame  rate 
from  a  single  facility  run  at  nominal  free  stream  Mach 
number  of  9.  The  field  of  view  is  again  ~  3.0  cm  x  6.5  cm 
and  the  nominal  free  stream  flow  velocity  is  ~  1,200  m/s, 
with  flow  direction  from  left  to  right.  The  nominal  static 
temperature  and  pressure  at  the  exit  of  the  fuel  injector 
is  ~  100  K  and  ~  140  kPa,  respectively,  based  on  a  CFD 
simulation  to  be  discussed  in  more  detail  below.  The  cal¬ 
culated  free  stream/jet  exit  momentum  ratio,  also  based  on 
CFD  simulation,  is  8.69.  The  corresponding  conditions 
of  the  free  stream  flow  directly  upstream  of  the  jet 
are  ~530  K  and  30  kPa,  respectively,  so  that  in  this 
instance,  the  jet  is  closer  to  pressure  matched  than  for  the 
images  shown  in  Fig.  7.  According  to  predictions  of  the 
CFD  simulation,  the  injected  jet  initially  expands,  with 
static  pressure  dropping  (due  to  mixing  with  the  free 
stream)  in  the  plume  to  ~  10-20  kPa,  followed  by  com¬ 
pression,  due  to  the  presence  of  a  rather  complex  structure 
of  oblique  shocks,  to  ~  50  kPa.  There  is  an  accompanying 
rise  in  static  temperature  to  300-700  K,  again  according  to 
CFD  simulation. 

The  downward  propagating  laser  sheet  was  brought 
incident  to  the  facility  from  top  right  to  bottom  left,  at  an 
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Fig.  8  10  KHz  NO  PLIF 
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angle  of  ~45°,  as  illustrated  in  Fig.  4.  The  gas  is  injected 
from  the  bottom-left  corner  at  30°  injection  angle  with 
respect  to  the  free  stream  flow  direction.  Note  that  the 
images  shown  in  Fig.  7  have  been  cropped  to  display 
only  ~400  (height)  x  900  (horizontal)  pixels.  The  aver¬ 
age  intensity  in  detector  “counts”  is  ~  1,000-3,000  out  of 
a  total  of  4,096,  which  was  achieved  using  an  intensifier 
gain  of  approximately  ~50  %  of  its  maximum.  The  cor¬ 
responding  background  noise  (mostly  detector  read  noise) 
is  ~20  counts. 

From  Fig.  8,  we  can  see  that  the  fuel  plume  injection 
rapidly  spreads  and  flaps,  filling  approximately  half  of  the 
total  30-mm  cavity  height  after  only  a  few  centimeters  of 
downstream  propagation.  In  actual  operation,  both  floor 
and  top  injector  banks  would  be  utilized,  which  would 
result  in  a  combustion  chamber  that  is  reasonably  mixed  at 
this  point.  It  can  be  seen,  however,  that  the  mixing  process 
is  highly  turbulent,  as  would  be  expected.  While  it  cannot 
be  determined  definitively  from  planar  imaging,  it  appears 
that  the  plume  spreads  rapidly,  and  turbulently,  in  three 
dimensions.  Of  course,  since  the  flow  speed  is  in  excess  of 
1  km/sec,  the  images  in  Fig.  8  are  temporally  uncorrelated 
at  the  10-kHz  frame  rate  in  which  they  were  acquired.  The 
uniformity  of  the  laser  sheet  is  also  a  potential  issue  when 
interpreting  these  images,  but  based  on  the  images  in 
Fig.  7,  and  recalling  that  the  laser  sheet  is  propagating 


downward  at  ~  45°,  it  appears  that,  with  the  exception  of 
the  edges,  sheet  non-uniformity  is  not  enormously  signifi¬ 
cant.  As  stated  previously,  only  the  central  two  inches  of  a 
nominally  four  inch  wide  sheet  was  used  for  the  actual 
imaging.  In  future  experiments,  it  may  be  possible  to 
measure  the  intensity  distribution  of  the  laser  sheet  on  a 
pulse-to-pulse  basis.  This  would  be  accomplished  by  use  of 
an  external  dye  cell  and  modification  of  the  imaging 
system. 

An  underexpanded  jet  structure,  similar  to  but  much 
smaller  than  that  shown  in  Fig.  7,  was  expected  to  be 
observed  in  Fig.  8,  near  the  bottom-left  side  of  each  image. 
No  such  structure  was  observed.  Three  factors  may  con¬ 
tribute  to  prevent  this  structure  from  being  observed.  First, 
the  laser  intensity  was  likely  considerably  lower  on  the 
edges  of  the  images  (it  also  varies  pulse-to-pulse, 
explaining  some  of  the  image-to-image  and  pulse  burst-to- 
pulse  burst  variation).  Second,  high  NO  concentrations 
were  probably  absorbing  the  laser  sheet  energy,  preventing 
fluorescence  near  the  lower  wall.  Finally,  it  is  possible  that 
beam  steering  due  to  mechanical  motion  of  the  facility 
resulted  in  some  clipping  of  the  edges  of  the  sheet, 
although  it  is  noted  that  the  beam  was  expanded  consid¬ 
erably  in  an  attempt  to  mitigate  this.  It  is  also  noted  that 
based  on  an  approximate  scaling  (using  a  CFD  simulation 
of  the  flow  field  to  be  discussed  below)  of  measured 
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absorption  in  a  room  temperature  gas  cell  filled  with  a 
mixture  of  5  %  NO  in  N2,  we  estimate  that  the  laser  beam 
absorption  in  the  current  experiment  could  be  as  high  as 
25  %/mm  in  the  central  portion  of  the  image,  corre¬ 
sponding  to  an  absorption  of  ~75  %  over  a  uniform 
~5  mm  “slab”  of  fluid.  Clearly,  this  is  a  high  level  of 
absorption,  and  future  NO  PLIF  measurements  should  be 
performed  with  lower  seeding  fraction.  Finally,  as  will  be 
described  in  more  detail  below,  quenching,  while  signifi¬ 
cant  in  the  jet  core,  is  not  sufficient  to  explain  the  lack  of 
observed  PLIF  signal  in  that  region. 

Examining  the  fuel  mixing  in  a  time-averaged  manner, 
the  upper  image  in  Fig.  9  is  the  average  of  40  individual 
NO  PLIF  images  obtained  over  several  facility  runs  at  the 
same  flow  condition.  The  nominal  per  pixel  spatial  scale  is 
again  ~75  microns/pixel.  From  this  image,  it  can  be  seen 
that  the  fuel  penetrates  the  primary  flow  at  the  injection 
angle  of  ~  30°,  propagates  at  30°  until  it  reaches  approx¬ 
imately  7  mm  in  height  above  the  floor,  and  then  turns  and 
propagates  approximately  parallel  to  the  primary  flow  axis. 
At  this  point,  the  plume  also  spreads  rapidly,  due  to  tur¬ 
bulent  diffusion.  After  approximately  2.5  inches  (6.4  cm) 
of  downstream  propagation,  the  plume  from  the  single  floor 
injector  has  penetrated  to  just  below  the  flow  centerline. 

A  parallel  computational  effort  was  performed  to  aid  in 
the  interpretation  and  analysis  of  the  experimental  data. 
The  bottom  image  in  Fig.  9  is  a  CFD  simulation  obtained 
for  the  same  flow  condition  as  the  NO  PLIF  average  image. 
The  simulation  was  performed  using  the  GASP  code  from 
Aerosoft,  Inc  (2000).  Although  the  actual  experiment  is 
highly  unsteady,  the  steady  state  Reynolds-averaged 
Navier-Stokes  (RANS)  simulation  performed  provides  an 
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Fig.  9  Comparison  of  averaged  NO  PLIF  image  with  CFD 
simulation 


“expected  mean  flow”  baseline  sufficient  to  aid  in  inter¬ 
pretation  of  data  without  the  computational  expense  of  a 
large-eddy  simulation.  The  main  features  of  jet  are  similar 
to  those  observed  by  Ben-Yakar  et  al.  (2006)  as  well  as  by 
Gruber  et  al.  (1995),  including  the  Mach  disk,  bow  shock, 
barrel  shock,  recirculation  zones  and  counter  rotating 
vortex  pair,  lending  confidence  to  the  accuracy  of  the 
overall  simulation. 

The  structured,  multi-block  grid  developed  with  the 
commercial  Gridgen  package  consists  of  approximately  48 
million  cells  and  exploits  symmetry  by  modeling  only  half 
of  one  injector  row.  For  the  simulation,  a  frozen-reaction 
chemistry  model  was  developed  from  existing  molecular 
thermodynamic  profiles  in  GASP’s  chemical  database  to 
predict  viscous  properties  of  the  helium  jet  and  nitrogen 
cross-flow.  Given  that  the  molar  fractions  of  the  NO  seed 
are  low  in  comparison  with  the  other  constituents,  the 
effects  of  the  NO  are  neglected  in  the  gas  model.  The  third- 
order  upwind-biased  Roe  scheme  is  employed  with  a 
weighted  essentially  non-oscillatory  limiter,  giving  an 
overall  third-order  accurate  solution  for  inviscid  fluxes.  To 
simulate  viscous  effects,  the  Wilcox  K-Omega  model  was 
used  with  a  compressibility  correction  to  account  for  over 
prediction  of  mixing  in  supersonic  shear  layers  (Aerosoft 
Inc  2000).  Turbulence  trips  were  modeled  using  a  forced 
intermittency  transition  along  the  walls  as  a  feature  of  the 
turbulence  model  implementation  in  the  GASP  code. 

To  simulate  the  effects  of  finite  sheet  thickness,  the 
computations  shown  in  Fig.  9  (bottom)  are  composed  of 
spatially  averaged  slices  in  the  spanwise  direction  (in/out 
of  plane  of  page).  Specifically,  the  image  represents  the 
unweighted  average  of  12  infinitesimally  thin  slices,  spaced 
52-microns  apart,  a  good  match  to  the  ~  600-micron  sheet 
thickness.  This  was  deemed  necessary  as  the  simulation 
indicated  a  highly  three-dimensional  flow,  with  significant 
streamwise  vorticity  in  the  vortex  pair  arising  downstream 
and  complex  accompanying  wave  and  shock  structure.  The 
intensity  scale  is  arbitrary  and  has  been  normalized  to  a 
peak  value  of  one  in  the  jet  core. 

For  comparison  purposes,  the  CFD  simulation  in  Fig.  9 
plots  the  relative  PLIF  signal,  5plif  as  follows: 

Splif  A  Nto,xNOfBQY  (1) 

where  Ntot  is  the  total  flow  number  density,  /NO  is  the  NO 
mole  fraction,  fB  is  the  Boltzmann  fraction  for  the  excited 
quantum  level,  and  QY  is  the  fluorescence  quantum  yield, 
given  by: 


where  A  =  5  x  106  s  1  is  the  spontaneous  radiative  decay 
rate  (Piper  and  Cowles  1986),  and  Q  is  the  collisional 
quenching  rate,  given  by: 
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Q(T)  =  21 A  x  I0~n NtotXm  (L)  1  (3) 


Note  that  Eq.  (3)  assumes  that  the  quenching  by  N2  and 
helium  is  negligible  compared  to  NO-NO  self-quenching, 
the  rate  for  which  at  300  K  is  taken  from  Paul  et  al.  (1996). 
This  is  a  reasonable  assumption  for  ~  95  %  of  the  mixing 
layer  given  the  5  %  NO  seeding  level  and  given  that  NO 
has  a  quenching  cross-section  that  is  almost  four  orders  of 
magnitude  higher  than  that  of  N2.  Eq.  (3)  also  assumes  a 
temperature-independent  quenching  cross-section,  as 
reported  in  reference  (Paul  et  al.  1 996),  which  results  in  the 
assumed  Tv 2  dependence  in  the  quenching  rate  due  to 
scaling  by  average  molecular  speed.  Note  that  the 
quenching  of  the  fluorescence  in  the  core  of  the  fuel  jet 
attenuates  the  fluorescence  signal  in  the  core  of  the  jet 
more  than  in  the  flow  in  the  duct.  Since  the  NO  injection 
pressure  is  very  high,  the  NO-NO  quenching  rate  is  also 
very  high,  ~8.1  x  108  s-1  at  the  ~140  kPa,  100  K  con¬ 
dition  at  the  injector  exit  where  yNO  ~0.05,  and 
4.8  x  107  s  1  at  center  of  the  field  of  view  of  the  jet, 
which  has  a  static  pressure  and  temperature  (according  to 
CFD)  of  ~37  kPa  and  500  K,  respectively,  and  yNO  has 
dropped  to  ~  0.025  due  to  mixing  with  the  free  stream 
fluid.  This  should  be  compared  to  the  NO  radiative  rate  of 
~5  x  106  s  1  (~200  ns  lifetime)  (Piper  and  Cowles 
1986).  Hence,  the  fluorescence  quantum  yield  varies  by  a 
factor  of  approximately  ten  over  the  imaged  region  of  the 
plume  and  is  10  x  smaller  in  the  plume  than  in  the  jet  core. 

While  it  is  difficult  to  make  an  absolute  comparison  of 
the  NO  PL1F  image,  which  has  not  been  corrected  for 
absorption  and  laser  sheet  non-uniformity,  and  the  CFD 
simulation,  it  is  clear  that  there  is  reasonable  overall 
qualitative  agreement.  One  significant  discrepancy,  noted 
above,  is  that  the  PLIF  images  do  not  show  any  indication 
of  the  barrel  shock  structure  predicted  by  the  CFD  in  the 
lower  left  region  of  the  jet.  As  described  above,  the  NO 
density  is  extremely  high  in  this  region  and  it  is  likely  that 
absorption  effects  cause  corresponding  reductions  in  signal 
that  mask  the  underexpanded  jet  structure.  In  future 
experiments,  lower  NO  concentrations  and  weaker  transi¬ 
tions  could  be  used  to  reduce  these  complicating  effects, 
thereby  producing  images  that  are  easier  to  interpret.  It  also 
appears  from  this  comparison  that  the  simulation  shows 
somewhat  less  spreading  with  distance  downstream  than 
the  experimental  data.  This  may  be  due  to  underprediction 
of  turbulent  mixing  as  well  as  the  absence  of  large-scale 
unsteady  structures  in  the  computation;  these  trends  are 
anticipated  with  the  RANS  approach.  Finally,  the  CFD 
simulation  predicts  a  thin  region  of  high  PLIF  signal 
starting  just  downstream  of  the  normal  shock,  where  the  jet 
fluid  turns  in  the  free  stream  direction,  as  well  as  a  more 
complex  structure  in  the  downstream  half  of  the  image. 


These  features  are  manifestations  of  the  rather  complex, 
and  three  dimensional,  structure  of  the  jet,  as  predicted  by 
CFD.  Both  of  these  predicted  features  are  not  observed  in 
the  experimental  image,  likely  due  to  turbulent  diffusion 
effects.  The  observed  maximum  intensity  in  the  center  of 
the  image  is  almost  certainly  due  to  sheet  non-uniformity, 
with  maximum  laser  energy  density  in  the  center  of  the 
sheet. 


4  Conclusions 

Nitric  oxide  PLIF  imaging  has  been  demonstrated  at 
10-kHz  repetition  rate  in  the  CUBRC  48-inch  Mach  9 
hypervelocity  shock  tunnel  using  a  pulse  burst  laser  sys¬ 
tem,  a  custom  Optical  Parametric  Oscillator,  and  an 
intensified  CMOS  camera.  During  the  run,  5  %  NO  in 
helium  is  injected  into  the  combustor  channel,  providing  a 
PLIF  fuel  tracer.  Up  to  ten  images  of  the  jet  mixing  in  the 
supersonic  flow  through  the  CUBRC  MURI  combustion 
duct  were  obtained  during  the  10-ms  test  time.  High  signal  - 
to-noise  images  were  obtained  under  typical  free  stream 
static  temperature  and  pressure  conditions  of  ~  530  K  and 
30  kPa,  respectively,  just  upstream  of  the  jet.  Comparison 
with  a  CFD  simulation  shows,  with  some  exceptions,  good 
overall  qualitative  agreement  of  the  predicted  and  observed 
penetration  depth  and  spreading  of  the  jet,  based  on  aver¬ 
aging  of  over  forty  individual  PLIF  images  obtained  during 
several  facility  runs.  This  work  demonstrates,  for  the  first 
time,  the  capability  of  increasing  the  data  rate  of  PLIF 
visualization  methods  in  high-enthalpy  impulse  facilities 
by  an  order  of  magnitude  utilizing  burst  mode  laser  tech¬ 
nology.  Future  work  will  emphasize  the  capture  of  more 
quantitative  flow  properties. 
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